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EQUATIONS FOR THE FIFTH SECANT VARIETY 
OF SEGRE PRODUCTS OF PROJECTIVE SPACES 


LUKE OEDING AND STEVEN V SAM 

Abstract. We describe a computational proof that the fifth secant variety of the Segre 
product of five copies of the projective line is a codimension 2 complete intersection of 
equations of degree 6 and 16. Our computations rely on pseudo-randomness, and numerical 
accuracy, so parts of our proof are only valid “with high probability”. 

1. Introduction 

Secant varieties have received growing attention in recent times, largely because of the 
fact that they provide a geometric model relevant to a wide variety of applications. The 
purpose of this note is to provide a case study in finding equations of secant varieties. For an 
introduction to secant varieties and their applications we invite the reader to consult [CGO] 
and the vast collection of references therein. The 5 factor binary secant variety is particularly 
interesting in light of recent results of Bocci and Chiantini [BC], that 2x2x2x2x2 
tensors are not identifiable in rank 5, but the generic tensor of that format has exactly 2 
decompositions. For > 6 factors, the binary Segre product is known to be ^-identifiable for 
most of the possible values of k below the generic rank, see [BCO] and [COV]. 

Fix a field k of characteristic 0. For % = 1,..., 5, let V) be a 2-dimensional vector space 
over k. Let V = &) 5 i=1 Vi. Let X be the 5th secant variety of nf= i P(b)) inside of P(V) (by 
the Segre embedding). The goal of this note is to prove * 1 the following statement. 

Theorem* 1 . 1 . The affine cone of X is a complete intersection of two equations: one of 
degree 6, and one of degree 16. 

We speculate that the homogeneous coordinate ring of any secant variety of any Segre 
product of projective spaces is Cohen-Macaulay. Theorem 1.1 confirms this for X. Using 
flattening and inheritance [L, Ch. 7], we get the following corollary: 

Corollary* 1 . 2 . Suppose n > 5. Let Vi,...,V n be vector spaces, let / 6 and /i 6 denote 
minimal generators o/crs^P 1 )* 5 ) and let F 6 and Fie, respectively, denote the natural liftings 
of fe and /j 6 to k[Vj <8) • • ■ <E> T4]. Let X be the fifth secant variety of 

P(W) x • • • x P(K) C P(Vi ® • • ■ ® V n ). 

Then the linear span of the GL(Li) x • • • x GL(T4) x Y> n -orbits of F 6 and Fie are equations 
that vanish on X. 

A flattening of a tensor A e V is a matrix constructed by viewing A as a linear mapping 
from the dual of one subset of the 5 vector spaces to the complementary subset. The basic 
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1 Our methods include probabilistic symbolic computations and numerical computations. Though they 
have been carefully tested and produce completely reproducible results, they are technically only true with 
high probability, or up to the numerical precision of the computers we use. To indicate reliance on such 
computations, we designate those theorems, corollaries, and propositions with a star. 
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fact is that if A has rank r, then a flattening F(A) has rank < r. So when non-trivial, the 
(r + 1) x (r + 1) minors of flattenings provide equations of secant varieties. In our case the 
only possible sizes of flattenings are (up to transpose) 2 x 16, or 4 x 8 , with maximum ranks 
2 and 4 respectively, so they do not provide non-trivial equations for tensors of rank 5. 

A next source for equations of secant varieties are exterior (or Koszul) flattenings (see, 
for instance, [CEO]) or, more generally, Young flattenings [LO] . For basic background, we 
invite the reader to consult [L], The basic idea of Young flattenings is to consider cases when 
the tensor product of a Schur module S M V := S w ki <g) S M2 V> 0 S /l3 Vj <g) S M 4 k 4 <S> S M3 V 5 with 
V contains another Schur module S^V, i.e., we have a linear map 

F /lM : S m V^S,V 

depending linearly on A e V. Subadditivity of matrix rank implies that if F^ I V {A) has rank 
p when A has rank 1 then F^^^A) has rank at most r ■ p when A has rank r. The art in this 
approach is to find good pairs of multi-partitions so that the dimensions of S M V and S^V are 
large with respect to p , so that the Young flattening has a chance to detect high rank tensors. 
In principle it is possible, but tedious, to list all possible Young flattenings and check which 
have a chance to provide meaningful equations. Because V) are all 2-dimensional, there are 
not too many choices for p and is, however our initial tries at finding Young flattenings that 
give non-trivial equations for the 5th secant variety were unsuccessful, and it seems that 
Young flattenings do not provide equations for this secant variety. 

After looking for Young flattenings unsuccessfully, we attempted a systematic search for 
equations via interpolation informed by representation theory. Our approach relied on com¬ 
puter calculations, which we explain in §2. The search for the equation of degree 6 is rigorous, 
but the search for the degree 16 equation is only correct up to high probability since we only 
show that it vanishes on sufficiently many pseudo-random points. Our search for equations 
was guided by our guess that this variety, having low codimension, would be defined by 
just a few equations, and because of the large symmetry group, that these equations would 
be semi-invariants. Our guesses are validated in §3, where we use these equations and one 
additional computer calculation on the degree of X (which is also only valid up to high prob¬ 
ability) to deduce Theorem 1.1. In §4 we provide a more detailed version of Corollary 1.2 
from a A-modulc (in the sense of [Sn]) point of view. 

In particular, the equations that we find provide modules of equations for all other 5th 
secant varieties of Segre products, both when the dimensions of the factors increase (by 
“inheritance”) and when the number of factors increase (by flattening). As far as we know 
these equations do not come from any known construction (such as Young flattenings), so 
they provide interesting new classes of equations for secant varieties. 

Acknowledgements. We thank Bernd Sturmfels for suggesting this problem, J.M. Lands- 
berg for helpful discussions, and Jon Hauenstein for providing Bert ini help. The software 
Bert ini [BHSW1], Macaulay2 [M2] and Maple were helpful for this work. Both authors ac¬ 
knowledge the hospitality of the Simons Institute for the Theory of Computing in Berkeley 
where this work was carried out. SS was supported by a Miller research fellowship. 

2. A SEARCH FOR EQUATIONS GUIDED BY SYMMETRY 

2.1. General idea. The variety X that we are studying has low codimension so we expect 
its ideal to be cut out by few equations. In addition, the defining ideal of X has a large 
symmetry group. If one polynomial is in the ideal of A', then so is the entire vector space 
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of polynomials in the span of its orbit. So we expect X to be cut out by a small number of 
semi-invariant polynomials. In this section we describe how we use all available symmetry 
to cut down our search for equations. 

Choose a basis e 0 , e 1 for V) so that we can identify the coordinates of P(V) with xi where 
I G {0, l} 5 . Let R denote the polynomial ring Sym(V) = k[xj | / G {0, l} 5 ] and let 
X := X(A) be the ideal of equations vanishing on X. Since R is graded, R = 0 d Rd, and we 
can compute X^ C Rd for each d. The most naive approach to determining Td is to evaluate 
a basis of Rd on dim Rd points of X using the parametrization of A", store the results in a 
matrix M and compute the kernel ker(M) = Td- In practice we use pseudo-random points 
on X with rational coefficients. 

In exact arithmetic, non-vanishing is a certainty, but vanishing might yield a false-positive. 
So this test gives an upper bound for dim I d and a probabilistic lower bound. The confidence 
in the lower bound may be increased by evaluating on more points of X. Alternatively, 
one can work over a function held over k of large enough transcendence degree (i.e., use 
parametrized points) where vanishing yields a genuine equation. The downside is that such 
computations are more expensive. This approach only works for small values of d as the 
dimension of Rd grows quickly. In particular, one can use it for d = 6 (the lowest possible 
degree in which the equations of cr 5 can occur, a basic fact from the theory of prolongation 
[LM, Corollary 3.4]) to discover the equation / 6 in §2.2, but it will not work for d = 16 (the 
largest degree we tested). 

However, R d has an action of the group SL^ 5 and we suspected that X is defined by 
invariants of SLA’. This gives a much smaller space of functions to search. Set Ud = 
(Sym (/ V ) sl 2 and Td = V® d . For each even degree d = 2 m the space of SL 2 "’-invariants in 
Td has a basis consisting of quintuples of Young tableaux each of shape (m, m) (there are 
no invariants in odd degree). A quintuple of tableau can be interpreted as a function on V 
by applying the Young symmetrizer associated to the quintuple of tableaux. Since Rd is a 
quotient of T dj it is spanned by linear combinations of quintuples of Young tableaux which 
now satisfy certain linear dependencies. 

An explicit basis of Ud C Rd may be found without explicit knowledge of all of the relations 
as follows. We can verify that a given set of quintuples of tableaux are linearly independent by 
evaluating them on dimension-many pseudo-random points of V. If the matrix constructed 
in this way has full rank, then we have a basis of that space of invariants. If not, we continue 
selecting random quintuples until a basis is found. 

Finally, the space Ud has an additional action of £ 5 . Assuming that there is a single 
minimal generator of X(A) in a given Ud, it must be a semi-invariant of £ 5 , so either an 
invariant or skew-invariant. Let Uf 5 (respectively Uf 5,sgn ) denote the subspace of Ud of 
£ 5 -invariants (respectively skew-invariants). In Figure 2.1 we list the dimensions of these 
spaces of invariants for degrees up to 16. The results follow from standard character theory 
calculations whose explanation we will omit. We focus our search for equations of X in the 
space of semi-invariants for SL(2) x5 x £ 5 . We describe this procedure in the case d = 16 in 
§2.3. The following proposition is a summary of what we found: 


Proposition* 2.1. There are minimal equations / 6 , / 16 vanishing on X of degrees 6 and 16. 
Both are invariant under SL^ 5 . Furthermore, fe is a skew-invariant under £5 while fi$ is 
a Tj^-invariant. 
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Degree d 

dim Ud 

dim Uj S 

dim f/f’ sgn 

2 

0 

0 

0 

4 

5 

1 

0 

6 

1 

0 

1 

8 

36 

4 

0 

10 

15 

0 

2 

12 

228 

12 

2 

14 

231 

2 

9 

16 

1313 

39 

10 


Figure 1. The dimensions of Ud, and its subspaces of S 5 -invariants and skew-invariants. 


To clarify: / 6 was constructed explicitly and we verified symbolically that it vanishes on 
X and that it is a skew-invariant. The polynomial / 16 was also constructed explicitly and 
verified to be a Es-invariant, but we only verified that it vanishes on a large collection of 
pseudo-random points on X, and hence it belongs to the ideal of X with high probability. 

2.2. The equation /g. Given a monomial in the Xi (the coordinates on P(V)), define its 
skew-symmetrization to be c -1 sgn(cr)a; (T (/) where c is the coefficient of Xj in the sum. 

The polynomial / 6 has 864 monomials and is the sum of the skew-symmetrizations of the 
following 15 monomials: 


— ^ 00000 ^ 01010 ^ 01101 ^ 10011 ^ 10100 ^ 11111 ) 
— ^ 00000 ^ 01100 ^ 01111 ^ 10011 ^ 10110 ^ 11001 ) 
— ^OOllO^OlOOO^OllOl^lOOOO^lOOll^lllll) 
^ 00100 ^ 01000 ^ 01111 ^ 10011 ^ 10110 ^ 11001 ) 
— ^ 00100 ^ 01010 ^ 01111 ^ 10001 ^ 10111 ^ 11000 ) 
— ^0010l3l010103l0111l31l0000 3110110^11001) 
— 3l001103l0100l3l0110031l000l31l001031lllll, 


3?000003?011003?0111l3?100103?1011l3?11001 > 
3-00000 3?011013?011103'100113?101103'11001 > 
3-00100 3l010103l011113^10000 31101113111001 > 
3l00U03l010003l0110l31l000l31l001031lllll, 
3l001003l010103l0111l31l001l31l010l31ll000) 
3l001003l0101l3l0111031l001l31l010l31ll000) 
3l001103l0100l3l0111l31l001l31l010031ll000) 
3-001113-010103-011013-100113-10100 3-11000 • 


There is an alternative description in terms of Young symmetrizers, following the same 
construction outlined in [BO]. The Young symmetrizer algorithm takes as input a set of fill¬ 
ings of five Young diagrams, performs a series of skew-symmetrizations and symmetrizations, 
and produces as output a polynomial in the associated Schur module. One can search over 
all possible shapes and fillings of tableaux for a fixed number of boxes and evaluate Young 
symmetrizers to find modules of equations in I(X). We describe this method more fully in 
the next section. In degree 6 the situation is particularly nice. It turns out that there are 5 
standard tableaux of shape (3, 3) and content {1,2,..., 6} and the following Schur module, 
which uses one of each of the 5 standard fillings, realizes the non-trivial copy of (S>]Li(^ 3,3^) 
inside of Sym 6 (V): 
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From this description of the invariant, one can use the classical symbolic method to verify 
that a general point of the fifth secant variety must be a zero of this invariant (see [R] for 
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complete descriptions of this type of argument, or also [O] ). This gives an unconditional 
proof that fa is in Z(A). Moreover, in §4 we explain how this description of fa also provides 
a generalization to 5th secant varieties of any larger number of Segre products of projective 
spaces of any dimensions. 

2.3. The equation fig. We follow the approach outlined in §2.1. Our search for new 
equations in degrees 8, 10, 12, and 14 did not yield any new equations, so we only describe our 
process in the degree 16 case. We used Maple, and the code for these computations may be 
found in the ancillary hies accompanying the arXiv version of this paper. Using the approach, 
we conclude that there are no E 5 skew-invariants in I 16 . We found that dim(Z D Uff) < 3 
and equals 3 with high probability. Note that for f G U w , we have / 6 /' G Uff if and only 
if f' G Uff” sgn since fa is a skew-invariant under S 5 . Since Uff’ sgu is 2-dimensional, we know 
that fa ■ U{f” ssn is a 2-dimensional subspace of Uff. Since dim(Z D = 3 (with high 
probability), we find that I has one additional minimal generator in degree 16. 

Let us describe in more detail the case of finding Eg-invariants, the skew-invariant case is 
similar. Let Qi, Q 2 , ... denote quintuples of Young tableaux all of shape (8, 8). The sum 

Fi := ^2 a -Qi 
o-es 5 

is in Uff . The Young symmetrizer algorithm (see [BO]) can be used to evaluate Fi on a 
point of V to test if it is non-zero. Each evaluation took between 500 and 23000 seconds 
and up to approximately 10GB of RAM on our servers 2 . 

We continued this randomized search for non-zero elements of Uff until 39 linearly inde¬ 
pendent invariants were found (the list of quintuples of fillings, data points, and code can 
be found in the ancillary hies). To verify independence, we chose 39 pseudo-random points 
Vi G V and verified that the matrix M := (F^Vj)) has full rank. While each evaluation 
Fi(vj ) is expensive, they are all independent computations. The distributed computation 
took approximately 2-3 days of computational time to verify the independence of F t . 

We then took 45 pseudo-random points pi G X and computed the matrix (Ff pj)) (only 
39 points are strictly necessary, but we included 6 more to increase the probability that our 
result is correct). This computation took an additional 2-3 days to complete. Finally we 
found that the matrix (F t (pj)) has rank 36, so there is a 3-dimensional space of invariants 
vanishing on (Py}y =1 an d vanishing on all of X with high probability. 

An expression of a basis of the kernel of the transpose of {F t {p 3 )) gives, in turn, a basis of 
Zi6 D Uif via linear combinations of quintuples of symmetrized sums of tableaux. Modding 
out by the space of invariants generated by / 6 , (in principle) one finds /i 6 . 

3. Proof of Theorem 1.1 
Lemma 3 . 1 . X has codimension 2. 

Proof. This follows from [CGG, Theorem 4.1]. □ 

Lemma* 3 . 2 . deg (A") > 96. 

2 One with 24 cores 2.8 GHz Intel Xeon processors with 144 GB RAM and another with 40 cores 2.8 GHz 
Intel Xeon processors with 256GB of RAM. 
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Proof. Computing the degree of an algebraic variety is a basic function of the software 
Bertini [BHSW1], see the ancillary hies and [BHSW2], We thank Jon Hauenstein for his 
help with this computation. The idea is to intersect X with a randomly chosen P 2 and 
find 96 points. One of the basic methods of numerical algebraic geometry (and Bertini, in 
particular) is numerical homotopy continuation, which relies on high precision numerical path 
tracking to follow the paths traced out by a homotopy from a system of polynomial equations 
with known roots to the desired system. For more details on these types of computations 
applied to tensor decomposition, see [DHO, HOOS]. 

The points found are represented by floating point numbers, so only satisfy the equations 
approximately. Though these methods have been rigorously tested through countless exam¬ 
ples, are open source and repeatable, there is still a chance that the computations yield a 
false positive result, so we only claim that the result holds with high probability, so the proof 
of the lemma may be read as evidence for the statement. □ 

Proof of Theorem 1.1. Using the notation from Proposition 2.1, let Y = U(/ 6 , fie), which is 
a complete intersection. We have ICf and deg A' > degU = 96. Since X is irreducible of 
codimension 2 (Lemma 3.1), and Y is equidimensional, Y is also irreducible (otherwise the 
degree inequality would be violated). So A' is the reduced subscheme of Y. Also, this implies 
that they have the same degree, so Y is generically reduced. Since Y is Cohen-Macaulay, 
generically reduced is equivalent to reduced. Hence X — Y is a complete intersection. □ 

4. Generalizations via flattening and inheritance 

As mentioned in the introduction, the equations / 6 and /i 6 provide modules of equations 
for all other 5th secant varieties of Segre products, both when the dimensions of the factors 
increase (by “inheritance”) and when the number of factors increase (by flattening). Inheri¬ 
tance for secant varieties was introduced in [LM], and was reinvestigated many times since. 
Inheritance can be encoded in the formality of A-modules [Sn] ; we hope that our description 
here provides an entry point to this formalism for the unfamiliar reader. 

The assignment (Vi,..., V 5 ) t—)■ Sym(Vi <g) • • ■ <g> V 5 ) is a multivariate polynomial functor 
R on 5-tuples of vector spaces, as is the assignment of the coordinate ring of the 5th secant 
variety of the corresponding Segre product. Similarly, there is a functor T t j defined by 

T id ( V u ...,V 5 )= Tor l Sym(yi 0 -® y 5 ) (i?(Ui, ..., V 5 ), k)j 

(the Tor groups are Z-graded, and the subscript denotes the jth homogeneous piece). When 
i — 1 , this is the space of minimal generators in degree j of the ideal of i?,(Vj...., V 5 ). 

As a representation of GL 2 , a SL 2 -invariant of degree 2 n is the Schur functor S nn . So 
Proposition 2.1 can be interpreted as saying that 83,3 E • • ■ E 83,3 appears with multiplicity 
1 in Tj 6 (coming from / 6 ) and S 8]8 Kl • • • E S 8]8 appears with multiplicity 1 in T 1)16 (coming 
from fio). The Koszul relation amongst fe and fie also shows that Snpi IE • ■ • E Snpi 
appears with multiplicity 1 in T 2 i 22 - Furthermore, if any other product of Schur functors 
S A i IE • • • IE S A 5 appears in any T, j then £(X k ) > 2 for some k since it vanishes when we 
evaluate on (k 2 ,..., k 2 ). 

There is another interpretation of our results using A-modulcs (see [Sn]): the 5th secant 
variety of the Segre product of projective spaces is a A-variety, and hence the assignment 
of a tuple of vector spaces to the space of degree cl equations vanishing on the 5th secant 
variety is a finitely generated A-module (which implies finitely presented using [SS, Theorem 
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9.1.3]). For d — 6, we have shown that S 33 IE • • • IEI S 33 (5 copies) are minimal generators of 
this A-module, and similarly for d — 16 and S 8 8 E • • • E S 8j8 . 
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